Dynamical fluctuations for periodically driven 
diffusions 



Navinder Singh 

Department of Applied Mathematics, Faculty of Sciences, H.I.T.- Holon Institute of 
Technology, Holon, Israel 

E-mail: navinder.phy@gmail.com 
Bram Wynants 

Institute of Theoretical Physics, K. U. Lcuven, Leuven, Belgium 
E-mail: Bram . WynantsOf ys . kuleuven . be 

Abstract. We study dynamical fluctuations in overdamped diffusion processes driven 
by time periodic forces. This is done by studying fluctuation functionals (rate functions 
from large deviation theory) , of fluctuations around the non-equilibrium steady regime. 
We identify a concept called traffic. This traffic, which was introduced in the context 
of non-equilibrium steady state statistics, is extended here for time-dependent but 
periodic forces. We discuss the fluctuation functionals of occupations and currents, 
and work out some specific examples. The connection between these and non- 
equilibrium thermodynamic potentials, their corresponding variational principles and 
their Legendre transforms, are also discussed. 



PACS numbers: 05.40.-a, 02.50.Ey, 64.60.-i,89.75.-k 



Dynamical fluctuations for periodically driven diffusions 



2 



1. Introduction 

Fluctuation theory is a cornerstone of statistical mechanics. On the one hand fluctu- 
ation theory naturally provides variational characterizations of the stationary regime. 
On the other hand, the fluctuations themselves are very interesting, especially because 
in recent years all kinds of mesoscopic systems have become observable and manageable. 
Think of e.g. molecular motors and transport through nanotubes in biophysics. Fur- 
thermore, and because of this, fluctuation theory has already helped in understanding 
nonequilibrium statistical mechanics, see e.g. jl] - [H]. 

In fluctuation theory one one has to distinct between static and dynamical fluctuations 
[13]. In static fluctuation theory one conditions on having a stationary regime in the 
infinite past, and measures the system at the present (time zero). This has been exten- 
sively studied in [31 0] for the hydrodynamic limit of several driven lattice gases. There, 
a Hamilton- Jacobi equation is derived for the fluctuations in the spirit of [2]. 

We will concentrate on dynamical fluctuation theory, as started by Onsager and Machlup 
in 1953 P^. In this theory, the system is continuously measured, and the fluctuation is a 
sequence of correlated non-equilibrium states (the time interval between two successive 
measurements is less than the correlation time). If the measurement time is infinite, the 
measured value of the observable will coincide with the 'typical' value. For finite times, 
averaged observables may or may not coincide with this typical value. If measurement 
times are large, then the fluctuations have to persist for a long time to be significant: 
one enters in the regime of large deviations. 

Extensions of the Onsager-Machlup theory have been made in recent years to 
nonequilibrium dynamics. We will continue in the lines of (TTJ [121 D3] • More precisely, 
this work is an extension of our previous work in [TTj . In [11] we have considered 
dynamical fluctuations of driven (by a time independent drive) overdamped diffusions. 
The most important results there are: 

1. The introduction of the concept of traffic, which complements the concept of entropy. 

2. The computation of the joint fluctuations of occupations and currents. 

3. A formulation of dynamical fluctuations in terms of thermodynamical potentials. 

4. An explanation of entropy production principles close to equilibrium. 

In the present contribution we carry on the program of [11] to time-dependent (but 
periodic) dynamics like in [15] . We will start here by introducing the model we work 
with. 

1.1. The model 

To introduce some notation, let us start with a simple model for an overdamped diffusion 
of a Brownian particle in d dimensions. Here and in the rest of this text we will always 
assume the process to be ergodic. It is described by the Langevin equation (using the 
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Ito interpretation): 



dx 



t — 



x(x t )VU(x t )dt + VD{x t )dt + ^2D(x t )dB t 



(1) 



Here, ci-dimensional vector describing the position of the particle. The d- 

dimensional space that the particle can move in is denoted by Q. B t is a <i-dimensional 
vector of standard Gaussian white noises, x is t ne mobility matrix and D is the diffusion 
matrix. Using the Einstein relation we have \ = (3D, where (3 is the inverse temper- 
ature of the surrounding fluid (bath). The potential U can be seen as the energy of 
the system (particle). We will assume that Q is either compact with periodic boundary 
conditions, or without boundaries. Think of e.g. a <i-dimensional torus or M. d . In the 
latter case the potential U should be sufficiently confining, which makes sure that the 
particle does not 'escape to infinity' (which happens e.g. for pure diffusion with U = 
in M. d ). An example of such a confining potential is a harmonic potential U — | Ylf=i x 1- 

The Langevin equation defines trajectories of the particle, a trajectory being the position 
of the particle as a function of time. Instead of considering trajectories, one can also 
take a probability density [1q(x) as an initial condition, and see how this evolves in time. 
This evolution is given by the Fokker-Planck equation: 



solves (EJ) with left-hand side zero, and is thus a stationary distribution for this process. 
We assume that U and x are such that this stationary distribution is unique and that 
for any p,Q we have that fi t — > Po for t — > oo (for more explanation and mathematical 
rigor, see [16]). Indeed, we are dealing here with a (simple) equilibrium process. This 
process we will use as our basic reference process. 

We add to the reference model ([T|) a time- dependent force f t . We take the time- 
dependence to be periodic with period r: ft = ft+r- The Langevin equation of the 
new stochastic process is now given by: 



where we shortened notation by defining F t (x) = ft{x) — VU(x). Note that the force 
which the particle undergoes depends on time in two ways: first of all because the 
force depends on the position of the particle and the position depends on time, and 
secondly because the force changes in time according to some deterministic periodic 
protocol. Note also that the system is now driven from equilibrium, because of the 
time- dependence of the dynamics. Our analysis includes the case in which the force f t 
is nonconservative, meaning there is no potential from which it derives, and as a result 
currents are generated. 




(2) 



The Boltzmann distribution 



Po(x) = \e~^ 



dx t = x(x t )F t (x t )dt + WD{x t )dt + y/2D(x t )dB t 



(3) 
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The corresponding Fokker-Planck equation for distributions p t is now: 




(4) 



with the probability current J Mt given by 

J ih( x ) = x(x)F t (x)fx t (x) - D(x)Vm{x) 



(5) 



As the forcing is now time-dependent, we do not expect to find a stationary distribution. 
But because the forcing is periodic in time, what we expect is that given a long enough 
time, measures p t will 'relax' to an evolution that is also periodic in time. In other 
words we assume that there is a measure pt{x) that solves the Fokker-Planck equation 
(jlj) and satisfies pt = pt+r, and moreover p t — > Pt for t — > oo. We will call this a non- 
equilibrium oscillatory state (NOS). Note that such a NOS becomes a non-equilibrium 
steady state (NESS) when the period of the forcing goes to infinity and the dynamics 
become time-independent. 

1.2. Questions 

As said before, when one considers particle densities, rather than the trajectory of one 
particle, the dynamics of the system is described by a Fokker-Planck equation, which 
is basically a deterministic equation. The determinism arises because of the law of 
large numbers, as a particle density by its very definition involves many particles. But 
for mesoscopic systems this assumption is not always valid, and fluctuations from this 
average "Fokker-Planck behavior" become important. The first question that then arises 
is: 

1. What are the right observables? 

These observables will be the empirical occupations p n ^ and currents j n ,t> where n stands 
for the number of periods (T = nr) that have lapsed and t G [0, r] for a specific time in 
each period. Then e.g. p n j,(A) counts the number of periods in which the particle is in 
a subset A of the total space at this specific time t within each period. We will explain 
this more exactly in Section I3TT1 

What we expect is that p n ^ —> pt and —> J pt when we let n — > oo. But, as the process 
is stochastic, for large but finite numbers of periods, anything can happen. There is a 
probability that, even after a long time, the system has not yet reached the NOS, but 
p n j resembles some other distribution, p t say. It is here that we enter the domain of 
dynamical fluctuations and large deviations [18j [19]. The probability of this happening 
will be exponentially small in time: 



where I{pt) is called a rate function, or a fluctuation functional. In the same way we 
can define I(j t ) for the currents, and I(p t ,jt) for the joint fluctuations of occupations 
and currents. It is here that the following and main questions arise: 



P{p n ,t = Pt) oc e 
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2. Can we calculate the fluctuation Junctionals? 

3. What is their thermodynamical meaning? 

Concerning question 2: in Section l3~2l we explicitly calculate I(fi t ,jt), see ( 12TT) . From 
this I(fit) and I{jt) can in principle be calculated through variational calculus, as 
I(fit) = infj t / (fit,jt) an d I(jt) = inf Mt I{^ti jt)- Explicitly calculating these is hard, 
but we work out some specific examples in the ensuing chapters. 

The main purpose of this text is to discuss the thermodynamic meaning of these fluctu- 
ation functionals. This is done in several steps at several places. First in Section [2] we 
define and discuss the main thermodynamic ingredients used throughout our argument: 
entropy and traffic. As we will see, entropy is time-antisymmetric and the less-known 
traffic is its time-symmetric counterpart. The latter is a measure of the dynamical ac- 
tivity of the process. Then, in Section [31 we keep track of these ingredients throughout 
the calculations. Additional thermodynamic interpretation is then provided in Section 
[6l where we try to make a connection with the idea of thermodynamic potentials. 

Finally, we study the limits of slow dynamics and small fluctuations. With a slow 
dynamics we mean here that the time-derivative of the forcing is small. This will be 
discussed in Section [3 

2. Entropy and traffic 

Our setting is a driven diffusive system which is not in equilibrium. Constantly, there 
is an entropy flux between system (the Brownian particle) and environment (the heat 
bath). This entropy flux is of course just a heat flux. If the system has at a certain 
time a density /i f , then the heat flux into the environment is equal to minus the change 
in energy of the system plus the work done on the system: 



Indeed, the first term on the right-hand side is the work done by the nonconservative 
force and the second is the energy change. Note further that this flux is a flux per unit 
of time (i.e. it is a rate). 

It has been known for a while now p2] that the entropy production can be defined and 
computed by comparing the probabilities of trajectories to the probabilities of time- 
reversed trajectories. This way of looking at entropy production is very convenient in 
dynamical fluctuation theory because it makes clear why it plays such an important 
role in the fluctuation functionals. On the other hand, in this way it also becomes clear 
that entropy is not the only important quantity. Entropy has its counterpart, called 
traffic (see [HI E21 US]), that becomes very important when going out of equilibrium. 
Therefore we will define and compute these quantities in this section, before going to 
the actual computations of fluctuation functionals in the next section. To do this, we 
first need to know how to compute probabilities of trajectories. 




(6) 
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2.1. Trajectories 

In statistical physics it is important to calculate expectation values of observables. These 
expectation values are actually averages over all possible realizations of the stochastic 
process at hand. For our model such a realization is the successive positions the 
Brownian particle visits during an interval [0,T], this is called a trajectory or path, 
which we denote by oj = (x t )o<t<T- The probability density of such a path is denoted 
by P At0 (a;), where fiQ denotes the density from which the initial position xq is chosen. 
Expectation values of observables O are then computed/defined by integrating over all 
possible trajectories: 

(O(u))^ o := J dojP, (oj)0(oj) = J dP, {u)O{u) 

Due to the Gaussian nature of the noise in the probability density of observing a 
path oj is proportional to[TT]. 

P^ (u) ex /io(x ) exp J q dt[x t - X F t - VD]D- l [x t - X F t - V£>] J 

where we have, for the sake of clarity, notationally suppressed the dependance of all 
functions on Xt- Actually, for the analysis in the rest of this text, we will use only relative 
probability densities. In this section we discuss the relative probability density of our 
process with respect to the reference process given by [TJ Denoting the path-probability 
density of the reference by P°, we get 

dP^o , \ _. fJ>o(xo) -am ( 7 \ 
dPV U) ■ Po(x ) [) 
where the action A is given by 

A{oj) = -t £ dx t °ft + ^l Q T dt[f t + |v - 2(VU)] X ft (8) 



in which the circle o denotes a Stratonovich-type stochastic integral (see Appendix B ). 
This relative density, described by the action (|Sj) is the key quantity needed to 
compute the fluctuation functionals in our framework. We closely examine its physical 
interpretation by splitting the action A into its time-symmetric and time-antisymmetric 
parts. We define the time-reversal operator 9 as follows: it reverses the trajectories in 
time 9oj = (xT-t)o<t<T, and it also reverses the protocol of the forcing 9 ft = fr-t- The 
anti-symmetric and symmetric parts of the action are then defined as 

S(u) = 9A(oj) - A{oj) = 13 T dx t o f t (x t ) (9) 

Jo 

7» = 9A{uj) + A(cu) = | jT dt[f t (x t ) + |v - 2(VU(x t ))]x(x t )f t (x t ) 

Physically, S(oj) is the excess entropy flux into the environment during the trajectory oj, 
excess with respect to the equilibrium reference process. We will name T(oj) the traffic 
in accordance with [IT]. 
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2.2. Entropy 

It has been known for a while now, that entropy production is a measure of irreversibility. 
In [TT] it was argued that the total entropy change in the world as a consequence of the 
process is given by 

dP^(uj) 



Stot(u) = log 



dP^Oiu) 

where /jLt is the measure that arises when /xq is evolved through time via the Fokker- 
Planck equation (J4]). Indeed, for an equilibrium process, like the reference process (pQ), 
this will be zero: 

S ^ u) = hg dP0l^) = ° 

This means that we can write 

dP^{u) . dP^6(u) 



S tot (u)= log -^H- log 



dP P0 (cu) ° dP po 9{u) 

= log . , - log — — - + S(u) 
Po{xo) Po{xt) 

with S(oj) the antisymmetric part of the action as defined in (Q. As said before, this S is 
the excess entropy flux into the environment. It is the work done by the force ft times /3. 

Let us compute the expectation value of Stot- The average of a Stratonovich integral 
is given in Appendix B see (IB. 81) : 



{Stot) a ■= I dPp{u)Stot{u) = S R (/j, ) - S R {iMr) + (3 I dt I dxf t J, 







in 



where Sr(h) is the relative entropy of the measure \x with respect to the reference 
equilibrium measure: 

/Lt(xj 
dxfi(x) log —-- 
Po{x) 

Using the explicit form of the reference equilibrium measure po oc e _/3C/ , we can rewrite 
S^t i n a wa y that is independent of the reference: 

{Stot)^ ■= J dP^ujStotiu) = s(fPr) - s(ju ) + J dt J dx[f t J M - fi t U] 
where s(/x) is now the Shannon entropy associated to the distribution /j, 
s(/i) = — dxfx(x) logfx(x) 



Compare this to ([6]). We see that the total entropy production is the change in shannon 
entropy plus the entropy flux between system and environment. 
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Using the Fokker-Planck equation we can still rewrite the total average entropy 
production in its shortest form: 

(Stot), = [ dtotbh) (10) 



with 



2.3. Traffic and its relation to entropy 

There is still a part of the action we have not discussed yet, that is, the time-symmetric 
part, or the traffic, second equation in fl9]). In contrast to the entropy production, this 
traffic depends very much on the reference process we take, so one should always keep 
in mind that traffic is an excess w.r.t to this reference. The average of the traffic can 
be computed and will be useful: 

( r io= en) 



with 



T t (jMt) := ^ t [ft ~ 2(W) + |v] X /t 



A straightforward calculation then connects this average traffic to the total entropy 
production: 

rt(pit) = l<T t QH) ~ \4 t= \^t) (12) 

where the second term on the right-hand sign is the entropy production in the reference 
equilibrium dynamics, but computed for a (nonequilibrium) distribution fi t . 



3. Dynamical fluctuations 

In this section we want to examine the fluctuation functionals that govern the asymptotic 
probabilities of observing nonstationary occupations and currents. For this it is very 
important to correctly define the observables we want to consider. 



3.1. Defining the observables 

To correctly define our observables we assume that the time of measurement T is very 
long, and a multiple of the period r: T = nr. We then define the empirical occupation 
density fi n ,t( x ) as follows: 

- n— 1 

f*n,t{ X ) = ~ y^J( X t+kr - X) (13) 

This density 'counts' for every time < t < r the number of periods in which the 
particle is at position x at this specific time within each period. In this sense this 
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empirical density is very detailed. As a function of x it satisfies J dxfi nyt {x) = 1. The 
reason for defining this detailed density is that it has the property that for any function 

u) t (x t ) 

— J dtw t (x t ) = i J dt J dxw t (x)ii njt (x) (14) 

As a consequence of ergodicity, we have that ^ J Q T dtw t (x t ) —> J dxw t (x)p t (x) for 
T — 7- oo. This shows us that for T — > oo we have that /i n t — > p t almost surely. In 
other words: 

lim V(y, n! t = fit) =0, n t ± p t 

T— voo 

lim V(fi n ,t = Pt) = 1 

T— »oo 

However, if /x f ^ p t , one can compute 

I(Pt) ■= ~ ^ lo S (^(^n,t = ^)) (15) 

where /(/it) is called a rate function, or a fluctuation functional. If it is finite, then we 
can schematically write 

V(fi n , t = fi t ) cx e~ TI ^ (16) 

by which we see that the probabilities of fluctuations fi t are exponentially damped in 
time. If I(fit) = oo, then the probability of the fluctuation fi t is damped even stronger. 
We will only consider fluctuations for which the fluctuation functional is finite. We have 
to admit that the equality sign in the probabilities V{fi n ,t = fit) should not be taken 
too seriously, the precise mathematical formulation of this and of rate functions can be 
found in [IB]. A more intuitive approach to large deviations is given in|19l 120] . 



The empirical occupation density is manifestly time symmetric. Its time antisymmetric 
counterpart is given by 

jn, t { X ) dt dXt+kr o5(x t+ kr ~ x) (17) 

k=0 

This measures the current at a; at a specific time < t < r in each period, and averages 
over n periods. By this definition we see that for any function wt 

Y dx t o Wt(x t ) = i J dt J dxw t (x)j n>t {x) (18) 

By this and ergodicity we can also see, together with formula ( IB. 81) . that for T — > oo 
we get jn )t —> J pt ■ In the same sense as for occupations, one can define the rate function 
I(jt). 

A more central and explicit starting point however is the rate function for the joint 
fluctuations of fi t and j t , defined by 

V{fl n , t = flu Jn,t = Jt) oc e~ TI ^ (19) 
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An important note before we go on calculating rate functions is that I(fi t ,j t ) = +00 

dt 



whenever [i t 7^ jj t+T or ^ + V • jt 7^ 0, see Appendix A Therefore we do not consider 



these cases. 

3.2. Joint fluctuations of occupations and currents 

To compute we use a standard large deviation technique, sometimes referred 

to as Cramer tilting, see e.g. [T8J HH]. Here is the recipe: modify the driving of the 
original dynamics, changing f t into some gt and we take care that this new and modified 
dynamics is chosen so that \i t and j t become both typical in the sense that 

^ + V- Ji = 0, jt = x(9t-VU)iM t -DVi h (20) 
which explicitly defines g t . 

We will prove that 

I ^Jt) = ^pJ dtrl((H)~J dtr t 9 (fi t )+^^ dt J dx(g t -f t )- Jt (21) 

The first two terms constitute an excess traffic averaged over one period of time (r). 
Note that from now on, we will add superscripts to quantities defining in which dynamics 
they are computed, i.e. in the original dynamics with / or in the modified dynamics 
with g. The second term is the average work over one period of time, done by the extra 
force g t — f t that is added to the original dynamics. Note that with the relation between 
average traffic and entropy production ( TI2"]) . we can also write 

I{ jh ,j t ) = ±-j dialog— J dta?(fi t ) + ^ dt J dx(g t -f t )-j t (22) 

To prove (12T|) . we write the probability of the fluctuations as 

P(V>n,t = IH,jn,t = jt) = J dPl (u)S(fI n>t - IH)S(jn,t ~ jt) 

dpf 

dP°M^^)Ktl n ,t - IHWnt - Jt) 

We use the formulas from section 12.11 to compute 
dpf 

)og- n g(L,)=A»(u)-A'{u) 

= l[T B (u)-T f H + Sf( U )-S»(u)] 

The expressions for entropy production and traffic are given in ([9]) and, and with (1141) 
and ([18]) these can be written as 

dpf t r 

lo § TT^(^) = TT / \ T !^n,t) - Tt(Pn,t) + 0(ft - 9t)jn,t] 



dP9 y ' 2r J 

- TI(fj lnjU j n - l 
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If we substitute this in the probability of the fluctuations, we get 

P(fin,t = fhJn,t = jt) = e~ TI ^ [ dP° (cj)6(Lln,t - lh)8{j n ,t = Jt) 



The second factor on the right-hand side is proportional to one in the asymptotic limit 
of T — y oo. This is because the average is now computed in the modified dynamics, in 
which n t together with j t represent the NOS, i.e. the typical behavior. This proves that 
the fluctuation functional is given by (I2ip . 



Using the explicit expression for r t (fi t ) (fTTl) . the fluctuation functional can also be 
written as 

t&tJt) = dt J dx Ut ~ J^ t D)-% - J M J (23) 

Indeed, the fluctuation functional for the joint probability of occupations and currents 
is fully explicit. To compute the fluctuations of only the occupations one 'simply' needs 
to integrate out the currents. Schematically: 

TI(vt) _ p(., , _ ,,\ _ / M, p -TI(p t ,jt) 



However, because we are considering the asymptotic limit of T — > oo, we are left with 
a variational problem: 

I(ji t ) = MI{m,j t ) (24) 

Vice versa we find the current fluctuations by taking the infimum over all fit of the joint 
fluctuation functional. 

4. Occupation fluctuations 

Let us try to compute the infimum in ( 1241) : the Euler-Lagrange equations give 

Jt = J* - (HxVt = fJttxlft -VU- Wt] - DVfit (25) 

where Vt is a Lagrange multiplier, making sure that + V-jj = 0. Unfortunately, there 
is no simple general solution to this equation. On the other hand, this equation gives a 
nice physical interpretation: to find the fluctuation functional for the occupations, one 
has to modify the original dynamics by adding a potential V t . This V t has to be such 
that it makes /x t typical, i.e. 

^ = -V • (fitxlft -VU- WJ - DV/it) (26) 
The fluctuation functional is then given by ( 12T1) with g t — f t — VV^: 
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This is thus equal to the excess in traffic plus the work done by the extra potential. The 
excess here is an excess of the original w.r.t. the modified process in which \i t is typical. 
Equivalently: 



I(jh) = ^J dt f dx^(VV t )x(VV t ) 



A downside of the results above is that there is no general explicit expression for I(/J>t)- 
Therefore we will consider some (simple) examples where one can get explicit results. 

4.1. One dimension 

In the case of a one-dimensional space Q, it is possible to explicitly calculate the 
occupation fluctuation functional. This is because the relation -j^ = —j' t (where the 
prime denotes derivation with respect to x) can be rewritten as 

for some number a G £1. This immediately gives the solution to ( 125*1) : 



-dy + j t (a) 



-J a 



+ /,(*) - U\x) 



dt 

Finally one can then determine jt{a) from boundary conditions. For example, let us take 
Q — R, and take a = — 00. if the potential U and the force are sufficiently confining, 
then jt(a) will be zero. This is for example the case when ft{x) — U'(x) is a polynomial 
in x with highest degree term equal to cx k with c > and k even. 
Another example is to take Q to be the unit circle. We take a = 0, and find jt(0) by 
using the periodicity condition on a unit circle V(0) = V(l): 

j.(Q) = \ ( f 1 f,( x )dx+ f 1 dx - I"" dx' d ^ X '^ 

Sl(Ht(x) X (x))-Hx \Jo Jo Vt(x) X (x) J 



dt 



4-2. Harmonic potential 

In this example we work in M d . We take the potential to be harmonic: U = k ^ i=1 x\ , 
and both the force ft and the diffusion matrix D to be independent of x. For any fixed 
time, this gives a detailed balanced process. It is the time-dependence of f t that pulls 
the system out of equilibrium. The corresponding Fokker-Planck equation is given by 

^ + V J Mt = with J M( = x(ft ~ kx t )/j, t - DV/j, t 
An advantage of this model is that we can explicitly calculate the NOS: 

yt z 

where the function at is given by 



a, 



dse -(t-s)k Xxfs 
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This means that the NOS is a Gaussian with a mean that is oscillating. Even for this 
harmonic potential one can not explicitly calculate I(fit) for a general /i t . However, we 
can calculate it for a restricted class of distributions. Inspired by the NOS we take fit 
to be an oscillating Gaussian: 



with b t a periodic function in time. A straightforward calculation then gives the solution 



Note that the current that minimizes I(/i t ,j t ) for this class of \i t is equal to -^^t- 
5. Current and velocity fluctuations 

As for the occupation fluctuations, a general explicit solution for the current fluctuations 
is not possible, even not in the one- dimensional case. Therefore we will again examine 
some simple examples. Before doing that, however, we will explicitly check the 
fluctuation theorem j9l [10]. For this we use the expression ( l23l) 



In the last equality, the second term within the square brackets integrates to zero, 
because the integration is over exactly one period of [i t . We are thus left with 



Z 



of ([H: 

Wt = ft-kb t - X - 1 %- 



so that the occupation fluctuation functional becomes 






and because the right-hand side does not depend on /i t we also have that 




By definition of the fluctuation functional, this is equivalent to 



P(jn,t = ~jt) 




which is an instance of the fluctuation theorem. 
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5.1. Divergenceless force 

We will consider here a dynamics in which V • (%i*t) = 0. In this case we immediately 
see from (J4]) that p t = -pr is just the uniform distribution, regardless of the period of 
the force F t . Note that we need to have a space with a finite volume in this case. 
Therefore we take Q to be the unit torus in d dimensions, which has volume \Q\ = 1. 
The corresponding NOS-current is J pt = x^t, which is then of course divergenceless. 
Let us therefore compute I(jt) for divergenceless currents. This means that we must 
minimize the joint fluctuation functional I(fi t ,jt) over all (time-) constant distributions 
fit = [i- This variational problem is straightforwardly solved, and one finds that [i = 1 is 
the minimizer for any divergenceless jt- This results in a current fluctuation functional 
that is quadratic in the current: 

Hit) = T~ A?* ~ xFt]D~% - X F t ] 



It 



'o 



5.2. Velocity fluctuations for a harmonic potential 

Here we will revisit the example of a harmonic potential, that was also studied for 
occupation fluctuations. While it is still difficult to find an explicit expression for the 
current fluctuations, it is possible to find explicit expressions for velocity fluctuations. 
The velocity associated to a distribution \i t and a current j t is v t = j t / ' \i t . 

We saw in the harmonic example for occupation fluctuations, that the current that 
minimizes the joint fluctuation functional for gaussian distributions, can be written as 
VtUt, with v t a factor independent oix. This factor is exactly the velocity. In this example 
we will consider the reverse problem: what is the fluctuation functional for velocities 
that are independent of xl For this we will first rewrite the fluctuation functional by a 
change in variables. The new probability density V for distributions and velocities is 

V(fi t ,v t ) = fj, t V([J,t, HtVt) 

But because the fluctuation functionals are defined as the logarithm of the probability 
densities divided by the time T in the limit for large times (as in (fT5]l ). we have that 

l{jJL u v t ) = Ifa, n t v t ) 

To find I(vt), we minimize the joint fluctuation functional over all /if. A straightforward, 
but somewhat tedious calculation gives that the minimizing distribution is a Gaussian: 

■ T,i( x i- b t,i) 



1 _M 



* = z e 



with the function b t equal to 

b t = I v s ds + f dt D~ l v t - 0f t + 0k [ v s ds 
Jo Pkrj [ Jo 

The fluctuation functional for velocities then becomes: 

Hvt) = ^ Jjx^vt - (ft - kbtMx'^t - (ft - kb t )] 
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6. A notion of thermodynamic potentials and variational principles 

In equilibrium systems it is useful to consider thermodynamic potentials (like free 
energy), as they have a clear physical meaning and characterize equilibrium via 
variational principles. Moreover, to go from one potential to another, Legendre 
transforms are used. Fluctuation functionals also bring with them variational principles. 
First of all, and mainly, minimizing them characterizes stationarity. This is easily seen 
for example for the occupations: 

V(fi n>t = fit) < V(fi n , t = Pt) oc 1 I(fi t ) > I(p t ) = 

In the same way, the joint fluctuation functional characterizes the typical occupations 
and the typical currents. In ( 12TT) we gave an explicit expression for this functional. Let us 
analyze this from the viewpoint of thermodynamic potentials and Legendre transforms. 
First of all, the first two terms in (1211) are time-averages of the traffic, or equivalently 
the entropy production ffT2l) . Let us consider the time-averaged traffic: 



H(fjL t ,ht) = \ f dtr ht (fi t 
P Jo 



where h t is an arbitrary force, which can be for example f t (the original dynamics of 
the system), or g t (the force that makes fit and jt typical). This functional H can be 
seen as a potential for the currents in the sense that the functional derivative of it with 
respect to the force, gives 
SH 



J h 



where J^ t is the probability current in a dynamics with a force h t , see (j4j). It is then 
natural to examine the Legendre transform of H: 



J dt J dxh t j t - H(n u h t 
The Euler-Lagrange equations to find the supremum give then 



G(nt,jt) = sup 

ht 



so that h t = g t is exactly the force needed to make the current j t typical, together with 
P-t- 

G(nu3t) = J dt J dxg t j t - H (fi t , g t ) 

Vice versa, if we take the functional derivative of G, we get: 
5G 

9t 



5j t (x) 

where gt is again the force that makes fit and jt typical. So G is a potential for the 
forces, just like H was a potential for the currents, and by Legendre transforms we can 
switch between the two. 
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We can easily rewrite I(fi t ,jt) in terms of G and H: 

G(jH,jt) + HQjLt, f t ) - J dt J dxf t j t 



Note that H(p t , 0) = 0, so that ^G(p t ,jt) is exactly equal to the fluctuation functional 
in the case that / = 0. In other words: 





I hi^hjt) = Io(jk,jt) + 2 



#(A*t, ft) - dt dxftjt 



This is nice, because the left-hand side is a fluctuation functional for a nonequilibrium 
dynamics, while on the right-hand side, the first term is a fluctuation functional in an 
equilibrium dynamics. The rest is thus 'the correction to equilibrium.' 

7. Small and slow 

7. 1 . Small fluctuations 

In this section we consider the regime of small fluctuations. In this regime a quadratic 
approximation is made of the fluctuation functionals, meaning that we are essentially 
in the regime of Gaussian fluctuations. More precisely, we take p t = p*(l + ep>i,t) and 
jt = Jp t + e ji,t with e a small parameter, so that \i t and j t are close to the stationary 
distribution and current. Up to second order in e, the joint fluctuation functional then 
becomes: 

I(fi t ,j t ) = ^- J dt J dx[ji >t - JptfJ>i,t + PtDV nx^ptD)' 1 

■ [ji,t - JptfJ>i s t + PtDVfii, t ] 

In this formula, one can see that the currents ji it and occupations \i\^ t are coupled be- 
cause of two reasons: first of all, because of cross-terms in the fluctuation functional, 
secondly because of the relation = —V ■ j\f 



The first of these reasons disappears in the regime of small driving, meaning that we 
will replace ft by ef t . First of all this means that p t = po(l + epi,t) and J Pt will be of 
order e. The fluctuation functional then further simplifies to 

I (( J 't,jt) = ^; J dt J dx[j 1) t(p D)- 1 j 1 j + p Vpi )t DVpi,t] (27) 
Indeed in this approximation there are no more cross-terms of currents and occupations. 

7.2. Slow dynamics and a comparison to time-independent dynamics 

Let us examine what happens if we take a dynamics which changes slowly in time. 
To parameterize this 'slowness' we take a small number e > 0, and change the time- 
dependent force ft to f et . This means that the period of the dynamics also changes to 
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r/e. What changes in the fluctuation functionals? To see this, note that the fluctuation 
functionals are time integrals over a period r, schematically: 

I —— J dtX(ht,ht) 

where X is a functional, which depends on some functions, here denoted by h t , but 
which can be fi t , ft, etc, (see the expressions for the functionals in the previous section). 
Making the dynamics slow by inserting e, we get 

e r /e i r 

I = - / dtZ(het, eh et ) = - dtX(h t , eh t ) 

T Jo T Jo 

where in the last step we just rescaled the integration variable t. So effectively, making 
the dynamics slow is equivalent to putting an e in front of all time-derivatives in the 
fluctuation functionals. As an example, consider the occupation fluctuations where we 
now get 

/(/it) = ^.f dtrft ^) ~ h Jo dtTft ~™ ( ' fJ * ) + ^So dt J dXfIt ^ 
In the limit e — > 0, a time-independent dynamics is recovered. For example, in the last 
formula this means that the last term drops from the fluctuation functional. As a check, 
in this limit the results of [H] are recovered. 



7.3. Small, close to equilibrium and slow 

If we consider the situation of (J27j) . and take also a dynamics that is slow (with the same 
e), then the coupling between the currents and the occupations completely disappears, 
up to second order in e. As a consequence, we can write 

I(fJ>t,jt) = I(fJ. t ) + I(Jt) 



with 



and 



I(^t) = ^ J " dt J dxptVfiijDVu^t 



c 2 



I(3t) = ^J dt J dxj^tiptD)- 1 ^ 

Moreover, because of this uncoupling one can write down a minimum entropy production 
principle for the occupations and a maximum entropy production principle for the 
currents, just like in the case of time-independent dynamics [11]. On the other hand, 
here so many approximations have been made that one could question the relevance of 
this. 



Dynamical fluctuations for periodically driven diffusions 



18 



8. Summary of the results 

In the calculation of the joint rate functional or fluctuation functional we have made use 
of the technique called Cramer's tilting, in which one changes the forcing such that the 
deviated observables become 'typical.' The joint rate functional consists of excess traffic 
and work done by the excess forcing both averaged over one period of forcing, see ( 12TT) . 
For the calculation of single rate functions (i.e. only for currents or for occupations) in 
the asymptotic limit of large T, this problem can be put in the form of a variational 
problem. For some specific cases the occupation and current functionals are derived 
explicitly. For the current rate functional the fluctuation theorem is verified. 

Importantly, in analogy with equilibrium thermodynamic potentials (like free 
energy) the non-equilibrium rate functionals (for example occupation rate functional) 
give us variational principles that lead to stationarity or the typical behavior. 

The time averaged traffic can be seen as a potential for the currents, as the 
functional derivative of it with respect to the corresponding force gives us the probability 
current (as explained in section IV). The Legendre transform of the potential for currents 
(the time averaged traffic) is the potential for the forces, i.e., the functional derivative 
of it with respect to the current gives us force. Thus Legendre transforms can switch 
between the two. Importantly, the joint rate functional can also be written in terms of 
this Legendre transform pair. 

Finally, we have considered the case of small fluctuations (the regime of Gaussian 
fluctuations). The joint fluctuational functional can be expressed in terms of the small 
occupation and current deviations. Interestingly, their coupling disappears in the regime 
of a driving that is both small and slow (section [7.3[) . As a cross check, we approximated 
the rate functional to the case of slow driving (the forcing period was extended to r/e, 
where e — > + ), and it reduces to the case of time independent forces, as expected. 
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Appendix A. Restrictions on fluctuations 

Here we will show that I{^ t ,j t ) = +00 whenever /i t 7^ /j, t+T or ^ + V • jt ^ 0. 

First of all, by definition (|T3|) we see that the difference between /i 4 and fi t+T is quite 

small: 

fi njt+T (x) = /i n ,t( x ) + -[$( x t+nT ~ x ) ~ $( x t ~ %)] 
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which is true for any u = {x t )o<t<nT- This also means that for any smooth bounded 
function w t we have 

dt I dxw t (x)[ii n j +T (x) - lM n ,t( x )] = _ / dt[w t (x t+nT ) - w t (x t )} 
J n Jo 

—7-0 for T — )• oo 

So we see that for any path, fi njt — fi nt t+ T — > in the distributional sense. Secondly, for 
any smooth bounded function Wt we have that 

dt dxw t [— + V • j nt t] = / dx t o-Vw t (x t ) / dtw t (x t ) 

Jo J dt n J n J 

Wq(xq) - W T {X T ) 

n 

— y for T — > oo 

So we see that for any path, + j' n t -> in the distributional sense. 

Al this means that fluctuations that do not satisfy fi t = fit+ T and + V • j n ,t = will 

have a probability that does not even survive in the asymptotic limit. 



Appendix B. Stochastic Integrals 

For the readers that are unfamiliar with the concept of stochastic integrals we give a 
brief description here, together with some useful formulas. The formulas in this ap- 
pendix resemble very much the formulas in the appendix of ref [llj . However, the fact 
that we have now time dependent forces complicates things a little, though not much. 

During computations in this article, one frequently encounters integrals over specific 
trajectories (paths) uj = (xt)o<t<T- As for diffusions these paths are nowhere 
differentiable, integrals like J dx t w t (x t ) are not well defined if one interprets them in the 
Riemann way. Instead there is room for different interpretations, the most common of 
which are the Ito and Stratonovich interpretations. 

ltd integral. For the Ito interpretation, the integral domain [0, T] is split up in a 
set of discrete points — to < t\ < . . . < t n — T, with Atj = tj — tj-i, such that 
At = maXj Atj — > for n — > oo. It is important to note that for diffusions as used in 
this article we have that for Atj — > 0, 

{Ax 3 ) 2 ah = (z» - xlJ{x\ - x b t .J 2D ab At 3 (B.l) 

where x a is the a-th component of the vector x, and D a t> is the corresponding matrix 
element. The stochastic integral interpreted in the Ito way is then computed as 

\ w t (x t )dx t = \im J2(x t -x t )w t (x t ) (B.2) 

J0 J = l 

Note that the function wt is a general time-dependent function. We do however assume 
that the time dependence of the function is smooth. This has as a consequence that it 
does not matter if we take uit J _ 1 or uit j in the definition of the integral. On the other 



Dynamical fluctuations for periodically driven diffusions 



20 



hand it does matter that one takes the evaluation point Xt_ 1 in the definition. 

For the Ito integral one cannot use the normal rules of integration. Instead one can 
check from (IB. 1|) and (1B.2j) that for any function w, 



T pT pT 

Vw t (x t )-dx t = w T (x T )-w (x )- dtw t (x t )- DV-Vw t (x t ) dt(B.3) 
'o Jo Jo 

where a dot means a derivative of the function with respect to time, and D'V ■ Vwt is a 

shorthand notation for Y, a ,b D abS^- 

Stratonovich integral. The Stratonovich interpretation differs from the Ito 

interpretation only in the points of evaluation of the function Wt- In this case Wt is 

evaluated in the midpoints of the time intervals: 

I w t (x t )odx t = lim - a?t,-_ 1 )wt( - „ — ) (B.4) 

J0 3=1 

where the symbol o is commonly added as a notation to distinguish between Ito and 
Stratonovich interpretations. Our definition of time reversal is that Xt —> xx-t and 
w t —> Wr-t- With this definition, the Stratonovich integral is time-antisymmetric. Note 
also that 

/ Vw(xt) o dx t = w{xt) — w(x ) — / dtw t (x t ) (B-5) 
Jo Jo 

so that the Stratonovich interpretation allows us to use the normal rules of integration. 

Relation between ltd and Stratonovich. It is easily found from flB.2j) and (1B.4|) that 

/ w t (x t )ox t = / w t (x t )dx t + / DVw t (x t )dt (B.6) 
Jo Jo Jo 

Averages of stochastic integrals. We can compute the average value of an Ito integral 

by using that x t solves the Ito-stochastic equation (J3J): 

/ wt(xt)dxt) =11 w t (xt)[xt(xt)Ft(x t ) + VD(x t )]dt) 

Jo ' Mo Vo ' 

T 

w t (x t )V2DdB t 

1 MO 

Observe that the last Ito-integral has mean zero since the integrand (evaluated at each 
mesh point x t ._^) and the increment B t . — B t ._ x of the Brownian motion are mutually 
independent, and the latter has zero mean. Hence, the mean value of the Ito integral is 



j w t (x t )dx t ) =(/ w t (x t )[x{x t )F t (x t ) + VD(x t )]dt 

Jo ' Mo \Jo ' 

= J dt J dxfx t w t [ X F t + VD} (B.7) 

The second equality uses that fi t , defined by the Fokker-Planck equation (TjJ, is the 
evolved density at time t when starting from /i at time zero. 
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Using the relation between Ito and Stratonovich, we also get: 

w t (x t ) o dx t ) =( [ [w t XtFt + w t VD + DVw t ]dt 



Mo Vo ' Mo 

T 







dt / dxw t J^ t (B-8) 
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